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Abstract 

We study the ground state phase diagram of the pseudospin model 
introduced by Doniach to describe the essential physics of Kondo lat- 
tices. We use variational trial states which augment the usual mean 
field solution by incorporating various intersite correlations. A com- 
posite spin correlation describing the antiparallel alignment of fluctu- 
ating triplets is found to be particularly favourable for large Kondo 
couplings. With this trial state, the magnetic-^to-Kondo transition 
is suppressed and the strong coupling ground state is ordered with 
strongly reduced moments. The relevance of the findings is discussed. 
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1 The model 



The basic question in the physics of heavy fermion materials is whether 
collective spin compensation can be taking place in a periodic array of lo- 
calized moments immersed in a conduction electron sea [Q. In case yes, we 
may speak about a collective Kondo effect even though it remains unclear 
to which extent the formation of an overall lattice singlet ground state can 
be likened to the single-ion Kondo effect. The most intriguing possibility is 
that spin compensation may go a very long way before it is stopped by the 
ordering of the residual tiny moments . 

As far as spins are concerned, the Kondo effect is just a spin compensation 
phenomenon. The emergence of a non-analytic energy scale in the impurity 
problem is connected with the existence of a large number of arbitrarily low- 
lying electron-hole excitations. Though variational methods indicate the 
existence of a lattice-coherence-enhanced Kondo energy scale for the nearly 
integral valent (Kondo) regime of the Anderson lattice p, ^, as well as for 
the Kondo lattice |^, they do not provide a proof that the ground state 
energy of the Kondo lattice contains non-analytic terms. The interplay of 
spin and charge degrees of freedom in the Kondo lattice may still prove to be 
quite different from what one has found for the impurity problem. It seems 
desirable to separate, if possible, the spin compensation aspect from all the 
other complications of "true Kondo physics". 

It has become accepted || that the competition of spin compensation 
and magnetic ordering can be described with drastically simplified models 
which contain just the spin degrees of freedom. The simplest of these is the 
Kondo necklace model introduced by Doniach 0. In addition to the localized 
spins S of the /-electrons we introduce a set of pseudospins r's which stand 
for the spin degrees of freedom of the conduction electrons. The number 
of S-spins is chosen equal to the number of r spins which implies that the 
possibility of a full spin compensation is a feature of the system. While it can 
be argued |TD| that a singlet ground state may arise at any band filling, it 
is more straightforward to associate the model with the Kondo lattice with 
a half-filled conduction band so that the number of pseudospins is equal 
to the number of conduction electrons. This has the additional motivation 



that the corresponding Kondo lattice has an insulating ground state [|TT| ; the 



appearance of a charge gap is a justification for omitting the charge degrees 
of freedom. 
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The simplest mean field treatment of the original necklace model yielded 
the beautiful result of a ground state phase transition from a magnetically 
ordered to a fully spin-compensated state |0. We are by now fully aware that 
strictly one-dimensional models are bound to show much subtler behaviour 
[^] but extending the model to the physically more relevant higher dimensions 
(no longer "necklaces" in the geometrical sense) makes us to expect that the 
mean field results are roughly correct. 

Our intention here is to improve the mean field approximation by includ- 
ing short-range correlations. Such an approach is expected to give sensible 
improvements over single-site mean field results in three dimensions where 
the mean field phase diagram should be qualitatively correct. However, in- 
termediate steps of our calculation can be executed free of further approxi- 
mations in one dimension. Therefore, we make the algebra (sum over local 
configurations) for the one-dimensional case but bear in mind that the char- 
acter of the results is meant for three dimensions. (In one dimension, the 
better estimate of the ground state energy is still believable but one should 
not trust the characterization of the ground state.) 

We study 5 = 1/2 Kondo necklace models 

H = jj2S,-n + W^E(rrrf+i + rfrf+, + pr/r^J (1) 
1=1 i=i 

The necklace is closed with the periodic boundary condition L + 1 = 1. 

Seeking correspondence with the insulating state of the Kondo lattice 
model would lead us to choose an antiferromagnetic Kondo coupling J > 0, 
and an antiferromagnetic intersite pseudospin coupling > 0. Actually, in 
Doniach's ^ original pseudospin model, the latter term was chosen to be 
purely x-y-like (p = 0), which gives a good imitation of propagating degrees 
of freedom. One should remember, however, that the underlying fermionic 
Kondo lattice problem had spin-rotational symmetry, and this has been lost 
by postulating the x-y form of coupling. We note that isotropic spin models 
(or anisotropic models with isotropy as a special case) have been discussed 
in the literature with the purpose of modelling Kondo lattice physics [|12|, 0]. 



In any case, a more complete understanding of the spin system makes the 
study of an extended model desirable. With this motivation, the study of (|^) 
with arbitrary signs of J and W, and with a general anisotropy < p < cxo 
is indicated. For reasons of convenience, we concentrate on the case p = 
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but we wish to emphasize that the method used here is equally applicable 
for p 7^ 0. 

Hamiltonians of the form (0) can be defined for either l^l > |r|, or 
\t\ > \S\, or \t\ = \S\, corresponding to under screened, overscreened, or 
exactly screened Kondo lattices. Generalizing Doniach's work on the exactly 
screened S = 1/2 model, we have discussed the mean field ground states of 



the underscreened [|13| and overscreened Kondo necklace models earlier. 
Here our main interest lies in going beyond the single-site mean field descrip- 
tion, therefore we confine our attention to the simplest case = 1/2, and 
|r| = 1/2. The Hilbert space of a lattice site is spanned by the four local 
basis states l^^r^) 

|1) = 11/2,1/2) |2) =1-1/2,1/2) 
|3) = 11/2,-1/2) |4) =1-1/2,-1/2) (2) 

The number of cases is reduced if we set p = 0. It is wellknown that 
the ferromagnetic and antiferromagnetic S" = 1/2 x-y models (on bipartite 
lattices) are physically identical since they can be connected by a canonical 
transformation. A bipartite lattice can be divided into alternate sublattices A 
and B so that nearest-neighbour bonds always connect different sublattices. 
Then the transformation 

U, = Y[ exp i-tnr;) (3) 
jeB 

changes the sign of the x-y term 

Ui fw^E(^^m + ^rl^ = -wj:{rM^, + rfrf^^ (4) 

\ i=l / 1=1 

A similar statement holds for our model: a 7r-rotation about the spin-^ axes 
for both the S"-, and r-spins on sublattice B 

f/2= nexpH7r(5| + r;)] (5) 

changes the sign of the intersite r-coupling while leaving the Kondo term 
unchanged 

\ i=l i=l / 
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1=1 



n 



W 



L 

El 

i=l 



'^i ''1+1 



(6) 



The second term on the right-hand side can be transformed to a non- 



interacting spinless fermion model ||T5[: for half-filhng, the ground state 



energy is —W/n. This provides a useful comparison for energy estimates in 
the weak-J regime. 

The inclusion of p 7^ would, of course, make the cases of positive and 
negative W genuinely different. Our variational method is, in principle, just 
as applicable for p 7^ as for p = 0. Restricting our attention to the case 
p = 0, where it is sufficient to consider < 0, is motivated by formal 
convenience: for ferromagnetic intersite coupling, homogeneous trial states 
can be used. 



2 Variational method 

We wish to map out the ground state phase diagram of H. For this purpose, 
we introduce variational trial states. Since a site i has two spins, the "true" 
spin Si, and the pseudospin fi which are coupled by the Kondo term, we 
can speak of the internal structure of the site, which is described as a linear 
combination of the four possible states 

\(t))i = a,\l)i + a2\2)i + a3|3)i + a^\4)i (7) 

In the single-site mean field theory, the internal state of the site i is taken 
to be independent of the instantaneous state of any other site; e.g., a trans- 
lationally invariant state would be described by the product wave function 

L L 

1$) = n i*^)^ = n(«iii)^ + «2|2). + «3i3). + «4i4)i) (8) 

i=l i=l 

In a ground state with antiferromagnetic long-range order, can be made 
sublattice-dependent 

L/2 
i=l 
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L/2 

= n(«^'|l)2i-i + «^|2)2i-l + «^|3)2.-1 + at\A)2^-l) ■ 

1=1 

(af |l)2i + ai\2)2i + ai\3)2i + af |4)2i) (9) 

but otherwise, the simple product form is retained. 

We propose to improve the variational description by including nearest- 
neighbour intersite correlations, according to the recipe: if site i is in the 
state and site z + 1 in the state then the amplitude acquires 

the additional factor 0(^j, ^i+i). For a translationally invariant state, we can 
write 

|vl/)=nAml$) (10) 

i 

where the intersite correlator 

4 4 

-Pj.i+l = X! X! C>i^i^^i+l) (H) 

has been introduced. 

The general structure of the variational trial state is that it is created by 
an intersite projection operator P acting on a mean-field reference state |$). 
On-site correlations (such as local singlet formation) are included in |$), 
and intersite correlations are controlled by P. Both |$) and P can contain 
variational parameters. 

Analogous trial states could be written down for higher-dimensional lat- 
tices. Actually, these would be physically more acceptable: our variational 
method is a correlated mean Geld method. Thus the overall appearance of 
our results is what we would expect for three-dimensional systems. Subtle 
features particular to one dimension |^ are likely to be missed by the present 
treatment. However, the algebra needed for the variational method can be 
executed fully in one dimension while it would involve further approxima- 
tions in higher dimensions. Therefore we stick to the one-dimensional case 
anticipating that care has to be exercised in interpreting the findings. 



The method we use has been introduced for spin chains by Virosztek [|T| 
later it was applied to the one-dimensional Hubbard model by Penc and one 
of us HI? 



We will be working with homogeneous states for which 

T.iS! + rn = (12) 
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The motivation comes from considering the single-site mean-field solution 
[^] of the original necklace model W > 0, J > 0, p = 0. In the simplest 
approximation, the ground state is an x-y antiferromagnet if J < W, and a 
collection of independent singlets if J > W, with a continuous phase tran- 
sition at J = W. The finding of a ground state phase transition between a 
magnetic and a "Kondo-like" state has long been the source of inspiration 
for continuing research in the Kondo lattice physics. On the other hand, it 
has been a matter of debate whether this transition is an exact consequence 
of the model ([l|), or an artefact of the simple approximation. At the simplest 
level, one can point it out that the description of the non-ordered state as 
strictly singlet is just a zeroth-order approximation since for any finite W/ J, 
however small, the process 

|l/2, -l/2)i| - 1/2, 1/2)2 |l/2, l/2)i| - 1/2, -1/2)2 (13) 

will mix in local triplets. 

It is natural to expect that the description in terms of local singlets be- 
comes correct only in the limit J/W — > 00, and for finite J, there will be 
spin-spin correlations between the sites. It is an interesting question whether 
the mean-field transition survives the inclusion of such correlations. 

3 Variational trial states for the necklace model 

3.1 Homogeneous states 

In accordance with (p!2D, in (§) we choose 

|cn| = |tt4| = and \a2\ = jaal = /? (14) 

One still has to decide the relative phase factors. 

The ground state wave function can be chosen as real, which says all the 
parameters are real. The remaining task is to specify the relative signs of 
the a's. This we do by requiring that non-diagonal processes should give a 
negative contribution to the energy, whenever that is possible. The on-site 
Kondo spin-flip term involves the factor 020(3'., the term can be made negative 
if 

sg(a2a3) = -sg(J) (15) 



7 



As we are going to see, there are two kinds of contributions coming from the 
intersite r-spin-flip processes. One of them, which acts between sites with 
antiparallel S'-spins, was illustrated in (|I^). Here all four local states appear 
once, so the term comes with the factor a\a20ij,a/^. The corresponding energy 
term can be made negative if we prescribe 

sg(aia4) = sg(J) ■ sg(iy) (16) 

In the other kind of hopping process, the S'-spins are parallel 

|l/2, l/2)i|l/2, -1/2)2 ^ |l/2, -l/2)i|l/2, 1/2)2 (17) 

In the corresponding term, all variational parameters are raised to even pow- 
ers, so the sign of the contribution is the same as the sign of W . In particular, 
for > 0, a homogeneous Ansatz does not permit to gain energy from this 
kind of process. 

No difficulties arise if we stick to the case p = 0: In (|^) we have shown 
that ly > and < are equivalent, so we can choose < 0, and have all 
contributions negative. Still, it is interesting to remember that the equivalent 
solution of the > problem is a two-sublattice antiferromagnetic state 
which can be generated from the homogeneous Ansatz by acting on it with 
the transformation IJ2 given in (^) 

L 

[/2-n(«|l).-/?|2)^ + /5|3).-a|4),) 

i=\ 
L/2 

= n("|l)2.-l - /^I2)2.~l + m2^-l - «|4)2.-l) ■ 
i=l 

{a\l)2i + l3\2)2i-m2i-a\A)2i) (18) 

The two-sublattice form is just Doniach's Ansatz for the case W > Q. 

Henceforth we keep W < 0. Depending on the sign of J, there are two 
different classes of wave functions. To be specific, we choose J > 0. The 
Ansatz for the ground state reads 

1^) = n + /?|2)^ - ^3\?>)^ + a|4),) (19) 
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The structure of the correlator can be represented by the matrix 



O 



( 0(6,6) 0(6,6) 0(6,6) 0(6,6) \ 

0(6,6) 0(6,6) 0(6,6) 0(6,6) 

0(6,6) 0(6,6) 0(6,6) 0(6,6) 

V 0(6,6) 0(6,6) 0(6,6) 0(6,6)/ 



(20) 



Considering the obvious symmetry 

0(66) = 0(6, e.) 



(21) 



and further restrictions arising from |-| symmetry, still leaves us with more 
variational parameters than one could easily handle. One has to try to guess 
what are the relevant correlations. We return to this question later. Before 
doing that, we outline the general formalism. 



3.2 Transfer matrix formalism 

We have to calculate the ground state energy 

{^\H\^>) 



E 



(22) 



and minimize it with respect to the variational parameters. 

Let us keep the matrix O and the choice of the a's as yet unspecified and 
calculate the norm (\E'|\E'). Expanding l^') in the orthonormal basis formed 
as the direct product of the single-site bases (0) 

1^) = E E • • • E n K0(6, 6+i)) I6)il6)2 . . . I6)l (23) 

we find that the norm is given by a sum over configurations, which is similar 
to the partition function of a one-dimensional classical lattice model. It is 
advantageous to introduce the transfer matrix T as 



^(6,6+i) = l«eJO^(6,6+i) I 



(24) 
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whereupon the norm becomes 

(^1^) = EE---E^(^i'6m6>^3)...T(eL-i,ama,6) 

= Tr{f'^)^x^ (25) 

where xq is the largest eigenvalue of the transfer matrix. 

The next step is the calculation of expectation values. The knowledge of 
xq suffices to determine the densities of quantities which are diagonal in the 
representation (|^). The density of sites in the state is 

^. = «e^ (26) 

The combined density of nearest-neighbour pairs in the configuration |^2), 
and its reverse |^i)|^2) is 

where the symmetry of the matrix O was exploited. 

The straightforward analogy with classical statistical mechanics ceases 
when we go over to the calculation of off-diagonal quantities such as spin- 
fiip amplitudes. At this stage, it becomes apparent that we are dealing with 
a genuinely quantum-mechanical problem. We follow the method introduced 
by Virosztek [Jl6| . 



Let us illustrate the method on the example of the spin-fiip part of the 
Kondo term acting at site m 

\ I m m ' m ml / 

m—l m~2 

E--- E EE E •••En 4- no^te.em) 

''^^mO{im-l-,im)0{^rn-,im+l){im\S'^T^ + S.^T^\^^ 

■a^>^0{U-i,C)OiCU+i)- n (40'te,^m)) (28) 

i=m+l 



The expression is rather like the norm (^51) except that at site m, a distur- 
bance has occured which, via the correlators, infiuences the sites m — l and 



10 



m + 1. (^) is still the trace of a product of matrices; however, in contrast 
to (p5D , not all matrices are T. The sites m — 1 and m + 1 are connected by 
the matrix K± rather than by 



'-\L—m 



(29) 



Doing the sums over and we recall that the spin-flip terms connect 
the states |2) and |3). It can be read off from (|28| ) that 

K±{im-i,ira+i) = 2 1 a^^^^ J ftatta 

■0(e„^-l, 2)0(2, U+i)0{U-i, 3)0(3, e^+i) (30) 

The exponential dominance of the trace by factors of Xq allows to deduce 
{^\s+T- + s-T+\^) {Xo\k±\Xo) 



Xn 



(31) 



where |Xo) is the eigenvector satisfying T\Xo) = xo\Xo). Because of the 
symmetry of T, it is of the form 



l^o) 



2(1 + A;21 



f k\ 
1 
1 

V k J 



(32) 



The intersite spin-flip process between sites m and m + 1 exerts an influ- 
ence also on sites m — 1 and m + 2. This can be expressed by a matrix M± 
which is analogous to K± 



m + 2 I 



■(0(1, 4)0(2, 3) [0{U-i. 1)0(4, U+2)0{U-i. 2)0(3, ^^+2) + 

0{U-i. 4)0(1, U+2)0{U-i. 3)0(2, U+2)] 

+0^(1, 3)0(e^_i, 1)0(3, U+2)0{U-i. 3)0(1, U+2) + 

0^(2, 4)0(e„_i, 2)0(4, U+2)0{U-i, 4)0(2, ^^+2)) (33) 



which has to be replaced into 



(34) 
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Another useful quantity is the transverse spin polarization which turns out 
to be the order parameter of the ground state. The matrix belonging to S!^ 
is 

iO{U-i, 1)0(1, U+i)0{U-i, 2)0(2, U+i) + 
0{U-i, 3)0(3, U+i)0{U-i, 4)0(4, U+i)) (35) 

while the matrix belonging to is the very similar 

iO{U~i, 1)0(1, U+i)0{U-i, 3)0(3, U+i) + 
0{U-i, 2)0(2, U+i)0{U-i, 4)0(4, U+i)) (36) 

The sign difference in (|35|) and (^) shows that (as expected for an antifer- 
romagnetic Kondo interaction) the 5* and r polarizations point in opposite 
directions. 

3.3 The minimization procedure 

Now we take particular forms of the Ansatz and work out the consequences. 
In the process, we hope to learn which correlations are most relevant for 
getting a correct description in different regimes of the coupling constant 
J/W. 

3.3.1 Simple spin correlations 

Since the mean field solution gives a magnetic-to-nonmagnetic ground state 
phase transition, the first idea could be to incorporate short range spin-spin 
correlations. Though only the r spins are subject to intersite interactions, 
their ordering induces a similar ordering of the S'-spins. (Since the Kondo 
coupling is antiferromagnetic, the local S'-moment is antiparallel to the local 
r-moment.) Therefore we incorporate in the Ansatz both r-r and SS 
intersite correlations, controlled by the independent variational parameters 
Xs and A^. Each antiparallel r-r pair brings a factor A^, and each antiparallel 
S-S pair a factor Ag. The corresponding matrix of the correlation coefficients 
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0(^j,^j+i) can be written as 



/I As A. A5A. \ 

A5 1 AsA^ A^ 

Ar AsA^ 1 A5 

V A5A. A, A5 1 / 



(37) 



In principle, we could have included intersite S — t correlations as well. 
Since the motivation for these is not immediately clear, and the minimization 
difficult enough with the parameters we already have, we omit them. 

Oi which describes fluctuating magnetism, should work reasonably well 
in the weak-to-intermediate coupling regime where it means an improvement 
over the mean-field finding of static long-range order. 

The largest eigenvalue of the transfer matrix is 



Xq 



a2 - /52)2(i + A|A2)2 + 4a2/52(A2 + X2y 



(38) 



and k appearing in the corresponding eigenvector 

1 



k 



2al3{\j. + A2 



/32)2(l + A2A2)2 + 4a2/?2(A| + A2)^ 



These quantities enter the various terms of the ground state energy. 
|-J, symmetries, the z-z part of the Kondo term can be derived from ( 

(^l^^^^l^) ^ (a^-/32)(l + A2A2) 

(^1^) 4^(a2 _ ^2)2(1 + A|A2)2 + 4a2^2(A2 + a2)2 

For the spin-flip part, we use (|3TD to arrive at 

(vl/|5+r„ + S„r+|vl/) 2/32AI A2(«A; + /?)2 



(39) 

Using 
Wi) as 



(40) 



Finally, (Q) is used to derive 

W ('\l/lr+r~ ,1 + r"r+ , 1 1^^ 

x^I'm m+l ' 'm m+l\^/ 

2 (^1^) 



-W 



xlil + k^) 
a2/j2A4(i + A2)3(afc + /3)2 

xm + k^) 



(41) 



(42) 



13 



To get a feeling for the structure of the result, let us first study the case 
J = when only the term (^2]) remains. Maximum freedom for r-spin- 
fiip is obtained for a = P = 1/V2 which leads to an expression which is 
independent of Ag: — 2A!^/(1 + Ar)^. This has its minimum at At = The 
minimum energy —8W/27 amounts to 94% of the exact value —W/tt. 

For small J/W, an expansion in terms of the small quantities jS/a — 1, 
and t — \/2 yields 

^H)^-^W--J-—C (43) 
^ ^ 27 9 361^ ^ ^ 

which corresponds to A5 ~ 1 and A,- ~ V^{1 — 3J/1QW). 

For general J/W, minimization was carried out numerically. The result 
for the ground state energy is shown in Fig. 1. In the interval < J/W < 1, 
the ground state energy shows an improvement over the simple mean-field 
result, which is quite substantial in the small- J regime. However, as J/W 
is increased slightly beyond 1, the solution seems to gradually approach that 
obtained by the single-site mean field treatment. 

We confirmed the existence of a sharp phase transition by a semi-analytic 
argument by expanding the energy in terms of the small parameter a/f]. 
Strictly for a = 0, the intersite hopping contribution (^) vanishes, and 
the Kondo energy has its minimum value —3 J/4 for A^At = 1. In the 
neighbourhood of the transition we expect (and find) that 6 = A5 — 1/ At 
is also very small, so we can make an additional expansion in terms of 6. 
The leading terms of the energy can be written as 

{H) ^ J + [MJ/W, K) + MJ/W, At) ■ 5] (^^^ (44) 

/2 is negative but an energy lowering due to an infinitesimally small S is 
impossible if /i is a finite positive quantity. We found that the minimum of 
/i changes sign at J/W ~ 1.059. For larger J's, the total energy increment 
in non-negative which requires a = 0, i.e., we are back at the single-site 
mean field solution. 

Below the threshold value, the ground state is ordered, having non- 
vanishing expectation values of and 

^ = ^''^ 
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Our experience with the Ansatz specified by p7| ) can be summarized hke 
this: introducing independent t-t and S-S correlations leads to considerable 
improvement in the description of the ordered state. However, there is still 
a phase transition from an ordered to a non-ordered ground state at the 
threshold value J = 1.05914^ which is near to J = ontained in the ordinary 
mean field treatment (Fig. 2). For J > 1.059iy, the description reduces 
to that obtained from the product trial state (P), i.e., it gives an array of 
decoupled singlets. The description of the high-J regime can be improved 
by postulating less obvious kinds of intersite correlations. 



3.3.2 Composite spin correlations 

Simple spin correlations failed to provide an acceptable characterization of 
the ground state for J/W > 1. To understand the nature of this state, let us 
remember that in the large- J limit, the process (|T3p prevents the system from 
freezing into a collection of singlets: it will keep on creating antiparallel pairs 
of local triplets. However, these pairs should then dissolve into singlets again, 
otherwise a high-energy situation remains sustained. Thus there must be a 
tendency for antiparallel components of local triplets to remain at nearest 
neighbour distance, which we try to enforce by 0(1,4) = 0(4,1) = t]. If 
we wish, we can assist the pair creation of local triplets by enhancing the 
antiparallel correlations in nearby singlets via 0(2,3) = 0(3,2) = (. Thus 
we are led to consider 

\ V 

It should be emphasized that here we are considering the correlations of com- 
posite objects made up of S- and r-spins. The correlations we introduced 
do not factorize into independent S-S and r-r correlations. 



1 1 r/^^ 
1 C 1 
C 1 1 
111/ 



(47) 
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The largest eigenvalue of the transfer matrix ( p^ is now 



(48) 



which belongs to an eigenvector of the form (|32D with 



The part of the Kondo coupling is 



while the spin-flip part is found to be 



(49) 



(50) 



(51) 



The r-spin-flip energy is given by 

_ ^ (^Ir+W + r^r+^Jvl/) ^ 2c.^/3^(l + r/Q [a(l + r/)A: + /j(l + Q]^ 
2 (^1^) x|](l + A;2)2 

(52) 

The energy expression can be put together from (^), (0), and (|I^). It has 
to be minimized with respect to the three independent variational parameters 
a/ 7], and (. 

Optimization has to be done numerically and the results are shown in 
Fig. 3 for the ground state energy and in Fig. 4 for (the absolute value of) 
the order parameter 

/c^x^ / .^ aP[Pil + C) + ak{l + r])f 



The most relevant feature is the suppression of the phase transition: the 
ground state is ordered for all J. The ordered moment is of 0(1) for weak 
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coupling and gets gradually suppressed in the strong coupling regime. We 
find it intriguing that a hint of reduced moment magnetism emerges in a 
model which is thought to correspond to the single-channel Kondo lattice. 
Details of the behaviour can be discussed in limiting cases: 
For J = 0, only the x-y term remains. Due to the symmetrical role 
played by local singlets and triplets in enhancing the mobility of the r- 
spins the minimum corresponds to a = /3, and f] = (. The lowest value 
is ~ — 0.281VF which we find at ~ 1.512. Thus even in the limit which 
is the exact opposite of what the trial state is intended for, a significant 
improvement over the mean field solution (energy —0.25W for 7^ = ^ = 1, 
a = (3) is achieved. The ordered moment is (r^) ^ 0.452. 

For large J, the density of local triplets has to be small, meaning a <^ 1, 
while the few triplets that are left must occur in antiparallel pairs so rj must 
become large. In contrast, the correlations governed by ( become relatively 
unimportant; for the sake of the present argument, we set C = 1- Introducing 
the convenient paramctrization a = cos(/9, j3 = sirup, Aip = 7i/2 — ip is one 
of the small parameters; the other is 1/r]. Expanding in these, the structure 
of the energy expression suggests to look for the asymptotic solution in the 
form 

ciW 
and 

A^='-f (55) 

where Ci and C2 are to be determined from optimizing the leading contribu- 
tions to the energy 

3 ^ /Cl 3ct ,2 



V = ^ (54) 



yielding 



and 



Cl = — z — i^V 




C2 = \hhrci (58) 
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The leading terms of the energy become 



which is obtained for 



and 



^ ' 4 32 J 



?7 ^ 8 1 




(59) 

(60) 
(61) 



The asymptotic behaviour of the order parameter is found to be 

(V3-1) W' 



1 + 



J 



(62) 



Thus our variational method recovered the correct order of magnitude ~ 
—W'^l J of the ground state energy in the large- J limit: it is what we would 
expect from perturbation theory. The nature of the problem is similar to that 
of the large-[/ behaviour of the Hubbard model where nearest-neighbour 
holon-doublon correlations were found to be important ||17|| . 

In retrospect we can identify the reason why simple spin correlations 
are insufficient in the large-J regime. The correlation matrix ( p7| ) enforces 
0(1,4) = 0(2,3) while with (g^) we have found that 0(1, 4)/0(2, 3) has to 
become very large as J /W increases. 



4 Discussion and conclusion 

We were trying to achieve a more detailed understanding of the behaviour of 
the Kondo necklace model, with the eventual aim of finding results which may 
be relevant to the physics of heavy fermion systems. We were considering 
the necklace hamiltonian (^ with p = 0, i.e., the simplest form introduced 
by Doniach 

The single-site mean field solution of ([l|) indicates that with increasing 
W/ J, a magnetic-to-Kondo (ground state) phase transition is taking place. 
While this is a physically appealing result, the obvious shortcomings of the 
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characterization of the high-J state as an array of disconnected singlets may 
lead to worries that the phase transition is merely an artefact of the approx- 
imation. 

We improved the variational description of the ground state by allowing 
for the presence of a variety of nearest-neighbour correlations. We carried 
out the optimization in the one-dimensional case where the transfer matrix 
method |16] can be used to calculate the relevant expectation values. 

In Section 3.3.1 we used the simple spin correlations whose presence can 
be inferred from the ordering. These led to a better ground state energy in 
the small-to-intermediate J regime, and pushed the phase transition point 
slightly upwards (Figures 1 and 2). However, the high-J state remained as 
structureless as in the simplest mean field approximation. 

In Section 3.3.2 we learned that the physically interesting ones are the 
composite spin correlations which could be expressed as expectation values 
of products of four spin-operators, involving both kinds of spins (we could 
choose {Sf S~^iT~^i)) . The relevance of these could be guessed from per- 
turbation theory: they describe that for J ^ W , the ground state is almost 
singlet, with a sprinkling of a few nearest-neighbour antiparallel triplet pairs. 
The inclusion of these correlations suppresses the phase transition completely 
(Figures 3 and 4), the ground state remains ordered for any finite J. Further- 
more, the ground state energy, and the concentration of local triplets have 
the order of magnitude expected from perturbation theory. The tail regime 
of the order parameter (Fig. 4) is a tantalizing hint that small ordered mo- 
ments may be a part of the physics of the orbitally non-degenerate Kondo 
lattice. 

However, we have to be extremely cautious about the conclusions to be 
drawn from our results. After all, for the one-dimensional Kondo necklace, 
powerful techniques have provided a number of essentially exact statements, 
and these tend to be in disagreement with our findings. General arguments 
suggest [0 that for almost the entire range of J /W values, the excitation 
spectrum is gapped but at some small J /W , the possibility of a phase transi- 
tion can not be excluded. In fact, exact diagonalization studies [|r^] revealed 
the existence of a Kosterlitz-Thouless type transition from the gapped high- 
J state (with exponential decay of spin correlations) to a gapless low- J state 
(with algebraic decay of spin correlations). In any case, the ground state 
never has true long-range order. 

This point in itself should not be too worrying. Our correlated mean 
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field approach is supposed to work where long-range order is in principle 
possible, i.e., foremost in three dimensions. It is just technical convenience 
which made us to stick to one dimension but we could argue that the general 
appearance of our results is the same as what a much more cumbersome 
three-dimensional evaluation should give. This would still allow us to hope 
that the three-dimensional pseudospin model would have a ground state with 
small ordered moments. 

We have to be, however, aware of a subtler kind of difficulty as well. Pre- 
scribing an Ansatz means that the system is permitted to seek a low-energy 
state in a certain manner. This can lead to a good estimate for the ground 
state energy (as it undoubtedly does) but does not necessarily imply that 
this is the natural way how the energy gain in question arises. To cite an 
example, a similar study of the Hubbard model [1^ gave the correct order 
of magnitude —t^/U for the energy at f/ ^ t but ascribed it to a metallic 
ground state which is patently false. Hence we should be warned that finding 
the (functionally) correct asymptotic form (|59|) of the ground state energy 
does not prove that the result (^) about the long-range order is basically 
right. In fact, preliminary results ||l9l obtained from the Oguchi approxima- 
tion indicate that the order parameter (r^) vanishes above a critical J /W. 
It remains an outstanding question whether the three-dimensional Doniach 
pseudospin model can support reduced moment magnetism. 
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Figure captions 



Fig. 1. Ground state energy (in units of W) versus J/W for a trial state 
with intersite SS and r-r spin correlations. Thick line: correlated mean 
field solution (transition at J/W ~ 1.059), thin line: simple mean field 
(transition at J/W = 1). 

Fig. 2. The order parameter (r^), for the same case as in Fig. 2. 

Fig. 3. Ground state energy (in units of W) for a trial state with composite 
intersite correlations, belonging to the antiparallel alignment of fluctuating 
triplets. Thick line: correlated mean field, thin line: simple mean field. 

Fig. 4. Order parameter (t^) versus J/W for the same case as Fig. 3. The 
ground state phase transition found in the simple mean field solution (thin 
line) is suppressed, according to the correlated mean field approach (thick 
fine), the asymptotic behaviour is ~ W/J. 
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